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ABSTRACT 

We have used data from the Herschel- ATLAS at 250, 350 and 500 to determine the 
far-infrared (FIR) properties of 50 Broad Absorption Line Quasars (BAL QSOs). Our 
sample contains 49 high-ionization BAL QSOs (HiBALs) and 1 low-ionization BAL 
QSO (LoBAL) which are compared against a sample of 329 non-BAL QSOs. These 
samples are matched over the redshift range 1.5 ^ z < 2.3 and in absolute i-band 
magnitude over the range -28 M l -24. Of these, 3 BAL QSOs (HiBALs) and 
27 non-BAL QSOs are detected at the >5ct level. We calculate star-formation rates 
(SFR) for our individually detected HiBAL QSOs and the non-detected LoBAL QSO 
as well as average SFRs for the BAL and non-BAL QSO samples based on stacking the 
Herschel data. We find no difference between the HiBAL and non-BAL QSO samples 
in the FIR, even when separated based on differing BAL QSO classifications. Using 
Mrk 231 as a template, the weighted mean SFR is estimated to be « 240±21 M Q yr" 1 
for the full sample, although this figure should be treated as an upper limit if AGN- 
heated dust makes a contribution to the FIR emission. Despite tentative claims in 
the literature, we do not find a dependence of C IV equivalent width on FIR emission, 
suggesting that the strength of any outflow in these objects is not linked to their FIR 
output. These results strongly suggest that BAL QSOs (more specifically HiBALs) 
can be accommodated within a simple AGN unified scheme in which our line-of-sight 
to the nucleus intersects outflowing material. Models in which HiBALs are caught 
towards the end of a period of enhanced spheroid and black-hole growth, during which 
a wind terminates the star-formation activity, are not supported by the observed FIR 
properties. 

Key words: galaxies: active - galaxies: infrared - galaxies: photometry - 
quasars:absorption lines 
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1 INTRODUCTION 

Mergers and other interactions are thought to trigger AGN 
which can lead to outflows that affect the chemical make-up 
of the interstellar and inter galactic media, and limit star for- 
mation in the host galaxy ( Bower et al.ll2006l; ICroton et al.l 
120061 ; iDe Lucia fc Blaizotll2007l ; iBooth fc Schavell2009l) and 
possibly in the large-scale e nvironment ( Rawlings fc Jarvisl 



12004 iMcCarthv et al.ll201oh . The role of these outflows has 
become an important feature of contemporary models of the 
formation and evolution of galaxies over cosmic time, be- 
ing required i n theoretical attempts to understand galax y 
'downsizing' (|Cowie et al, 1 1 19961 ; IScannapieco et al.1 iBoEl). 
In se mi-analytic models of galaxy assembly ( Granato et al.l 
12004 ) AGN feedback removes dense gas from the centres of 
galaxies, heating up the surrounding Intergalactic Medium 
(IGM) , while at the same t ime enriching the I GM with met- 
als. Studies such as that of lGabel et al. |2006l ) indicate that 
outflows have high metallicities which could serve to enrich 
the IGM. Due to the IGM's low density, however, it cannot 
cool efficiently and so cannot fall back to the galaxy to fuel 

star formation. 

Simulations of galaxy mergers (|Springel et al.1 l2005al ) 
show that AGN feedback can heat the Interstellar Medium 
(ISM) and inhibits star formation, with outflows provid- 
ing the energy and momentum feedback for the ISM of 
the host galaxy. Efforts to include th ese mechanisms (e.g., 
ICroton et al1l2006l ; ISiiacki et al]|2007l ) have focussed on two 
modes: a 'radio mode' whereby a relativistic jet can heat in- 
terstellar gasj amijritraclusl^^ further infall 
of mass !|Bestll2007l ; iBest fc Heckman|l2012l ). and a 'quasar 
mode', which serves to stop star formation with an out- 
flow of greater mass but with a lower velocity driven by the 
radiation from the quasar. At the same time, the quasar 
mode removes any leftover gas that might serve to enshroud 
the galaxy, turning it into a 'classical' Quasi-Stellar Object 
(QSO). Recent work finds that when feedback is incorpo- 
rated, the ACDM model provides a good fit to observed 
galaxy s tellar-mass fractions and the _R"-band luminosity 
function (|Bower et al.ll200d ; ICir"asuolo et alll2010l ). However, 
if feedback is excluded in galaxy group simulations, the tem- 
perature profiles are highly peaked and in disagreement with 
many of the observed properties of galaxies. The simulations 
also suffer from the well known overcooling problem, the re- 
sulting st ellar mass fraction bein g several times larger than 
observed (|McCarthv et al.ll20ld ). 

DD 

Gas outflows from AGN are primarily detected in X- 
ray and ultraviolet absorption against the inner portions of 
the accretion disk and/or the more extended broad-line re- 
gion. Indicators for the origin of an absorption line system 
are (1) velocity width, (2) partial coverage, (3) time vari- 
ability, and (4) high metallicity. Historically, the criterion of 
velocity width has led to the use of three categories of intrin- 
sic absorbers. Those with the greatest velocity dispersions 
are termed 'Broad Absorption Line' Quasars, or BAL QSOs 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. The H- ATLAS website 
is http://www/h-atlas.org/ 

tE-mail:jclOacd@herts. ac.uk 



(e.g lWevmann et al.ll 199 J ) . These have absorption lines with 
FWHM greater than 2000 kms -1 . Those with very narrow 
absorption lines are kno wn as 'Narrow Absorptio n Line' 
QSOs, (NAL QSOs; e.g. lHamann fc Ferlandlll999l ). These 
have FWHM of less than 500 kms" 1 . Finally the last type 
are known as Mini-BALs. These are those intrinsic absorbers 
whos e velocity widths are between 500 and 2 000 kms -1 
(e.g., lHamann et al.lll997l . IChurchill et al.lll999l ). 

There is still some debate in the literature regard- 
ing BAL QSOs, with the observed fraction (~ 15 per 
cent) being attributed to orientation, an evolutionary phase 
or some mixture of both. In the orientation hypothe- 
sis, BAL outflows are present in all QSOs but they are 
only viewed as BAL QSOs when our line-of-sight inter- 
sects the solid angle subtended by the BAL outflow. In 
this situation, high column density a ccretion disk winds 
are accelerated by radiation pressu re (I Murray et al.l 1 19951 : 
iMurrav fc Chiangj[l99e1 : lElvisll2000l ). This model attributes 
the fraction of QSOs to the fractional solid angle coverage 
of the BAL regions. Furthermore, it is able to account for 
the similarity of the continua and emission line features in 
BAL QSO and non -BAL QSO spectra once reddening is 
taken into account (|Wevmann et al.1 Il991 : iReichard et al.l 
2003al) and fits well with unified models ( Antonuccil 1993k 



Elvis! |2000| ). Mid-infrared and far-infrared studies also 



show no difference in emission, consistent with there be- 
ing similar dust masses within the host galaxy and BAL 
QSOs being an orie ntation effect (jGallagher et al.l 120071 ; 
lLazarova et al.ll2012l). However, spectropolari metry studies 
jQgle et al.lll999l ; lDiPompeo et al.ll2010M201ll ) indicate that 



BAL QSOs are seen at a wide range of inclinations, which 
suggests that BAL QSOs may not be just a simple orienta- 
tion effect. 

The evolutionary scenario is rather different. In such a 
model, BAL QSOs are young QSOs that may have recently 
undergone a merger/starburst, and are observed as BAL 
QSOs when two special conditions are met; firstly, the QSO 
luminosity has reached a sufficient strength to accelerate gas 
to several thousand kms -1 and secondly, there is a large 
mass of diffuse gas and dust in the nuclear regions leading 
to a rapid mass-loss phase. These BAL outflows could also 
play a role in terminating star formation and removing much 
of the obscuring dust and gas le ading to a classical optically 
and s oft X-ray luminous QS O (|Voit et al.l 1 19931 ; lHall et al.l 
120021 ; IPage et al.ll2004 l201ll ). If BAL QSOs are caught to- 
wards the end of an epoch of enhanced star-formation activ- 
ity then high dust masses within the host galaxies of BAL 
QSOs will yield higher average FIR flux densities. Since far- 
infrared emission is optically thin, there should be no de- 
pendence on the source orientation. This difference in the 
FIR properties in these two hypotheses can in principle be 
used to distinguish between them. 

lOmont et all l| 19961 ) and ICarilli et all (|200lf ). both 
working at frequencies around 250 GHz, showed that there 
was weak evidence that BAL QSOs were more luminous 
at millimetre wavelengths than ordinary QSOs. This would 
favour an evolutionary scenario, but the work did not con- 
centrate exclusively on the BAL QSO phenomenon, or com- 
pare star formation rates, and the samples w ere small (4 
and 7 BAL QSOs respective ly). Other work by [Lewis et al] 
(|2003l ). IWillott et all (|2003l ) and iPriddev et al. I (|2007h was 
specifically aimed at improving our understanding of the 
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BAL phenomenon, with larger samples of 7, 30 (41 non-BAL 
QSOs for a compariso n sample) and 15 BA L QSOs (supple- 
mented by data from IPriddev et al.l 120031 ) respectively. All 
of these samples were observed using the SCUBA camera 
on the James Clerk Maxwell Tel escope at wavelengt hs of 
450 and 850 ym. Th e samples of IWillott et all (|2003l ) and 
IPriddev et ail l|2007h are composites of previous work by the 
authors at FIR wavelengths, and therefore have non-uniform 
sensitivity coverage. None of these studies find strong evi- 
dence for BAL QSOs to be more submillimetre luminous 
than non-BAL QSOs, a finding consistent with the orien- 
tation model for BAL activity. Within the B AL QSO pop- 
ulation itself, however, IPriddev et all |2007l ) find tentative 
evidence for a dependence of submillimetre emission on C IV 
absorption-line equivalent width and discuss the implica- 
tions within the framework of simple evolutionary models. 

In this paper we use Herschel jPilbratt et alj l201Ch 
data to study the FIR properties of matched samples of 
BAL and non-BAL QSOs allowing us to build on previous 
work with improved statistics, better wavelength coverage 
around the FIR dust peak and better matching of the con- 
trol sample. Throughout this work we use a cosmology where 
Ho = 70 kms -1 Mpc~\ fi m = 0.3, and Q A = 0.7. 



2 THE DEFINITION OF A BAL QSO 

BAL QSOs have broad absorption lines which arise from res- 
onance l ine absorption in gas outflowing with v e locities up 
to 0. 1 c JWevmann et al]|l99lt lArav et all l200ll ; lHall et all 
120021 ; iReichard et al.ll2003al V As a subclass of AGN, BALs 
until recently were thought to make up around 10 — 15 



per cent of the sources in QSO surveys dWeymann et all 



;yn 

199ll; iTolea et alj |2002| ; IReichard et~afl l2003ah . However, 
Trump et alj (|2006l )~ quoted a much higher BAL fraction, 



i.e., an observed fraction with in the 3rd Data Release of the 
SDSS (|Schneider et al.ll2005l ) of 26 per cent. It is possible 
that selection effects may bias us against BAL identifica- 
tions, mea ning that the int r insic fr action is far higher than 
obser ved jHewett fc Foltzl 120031 ; iDai. Shankar fc Sivakofll 
120081 ). Perhaps the most widely accepted classification 
scheme for whether a QS O is a BAL QSO is th e Balnic- 
ity Index or BI, defined by I Wevmann et alj (|l99ll ). 
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Here f(v) is the continuum-normalized spectral flux at a 
velocity v (in kms -1 ) from the Civ line rest wavelength 
1549 A (in the system frame). The dimensionless value C is 
unless the observed spectrum has fallen at least 10 per 
cent below the continuum for a velocity width of at least 
2000 kms -1 in the absorption trough, at which point C is 
set to 1. Traditional BALs are defined to have BI > 0. 

We construct our BAL sample using BAL QSOs se- 
lected using the BI metric of equ ation [T] and a versi on of 
the metric which is 'extended' as in I Gibson et ail l|2009h . the 
lower limit of the BI being kms -1 . 

BAL QSOs can be split into subclasses based on their 
observed spectral features. The first of these classifica- 
tions are High-ionization Broad Absorption Line (HiBAL) 
QSOs. These objects contain absorption in LyaA1216A, 
NvA1240A, SilvA1394A and CivA1549A. They are the 



most prevalent sub-population, making up around 85 per 
cent of BALs. The Low-ionization Broad Absorption Line 
(LoBAL) QSOs contain all of the previously mentioned 
absorption features seen in HiBALs, but also contain ab- 
sorption features in MgllA2799A and other low ioniza- 
tion spec ies. These objects compr ise around 15 per cent 
of BALs jSpravberrv fc Foltzlll992l estimated 17 per cent). 
HiBALs and LoB ALs are found to have r edder continua 
than n on-BALs (iBrotherton et al. I |2001| ; IReichard et al.l 
l2003bl ; IScaringi et al.1 l2009h . with LoBALs also being 
significantly redder than HiBALs jWevmann et al.1 Il99ll ; 
ISpravberrv fc Foltzl 1 19921 ). The third and final class are 
FeLoBALs, which show the absorption features of LoBALs 
as well as absorption features arising from metastable ex- 
cited levels of iron. The effects of these definitions on each 
measured variable are discussed in greater detail in Section 
4. 



3 DATA AND SAMPLE SELECTION 

The BAL QSO sample that we use in this paper is a small 
proportion of those objects identified as a BAL QSO within 
the SDSS QSO Catalogue (|Gibson et al.ll2009l ). which num- 
ber 5039 BAL QSOs. By using the SDSS catalogue (both for 
BAL and non-BAL samples) a redshift of z < 2.3 provides an 
upper redshift limit where the SDSS colour sel ection criteria 
are r elatively unaffected by BAL absorption (|Gibson et al.1 
I2009D . We use absolute magnitude in the SDSS i band since 
it is les s affected by dust re ddening than the bluer bands. 

As lGibson et""aH |2009l ) used the BI as a classification 
mechanism, it is possible that some BALs will have been 
missed from the sample, and others included within the cat- 
alogue will be mis-classifications (though the strictness of 
the BI should decrease the number). 

In this paper we use observations from th e ESA Her- 
schel Space Observatory (|Pilbratt et al.l |2010| ) as part of 
the Hersch el Astrophysical Terahertz Large-Area Survey 
(H-ATLAS, lEales et al.ll2010h survey phase 1 data set, re- 
duced in the same w ay as the SDP d a ta set and with the 
same 5 a flux lim its (jlbar et al.1 l20ld ; IPascale et all |201ll : 
iRigbv et al]|201ll ). It is the largest Herschel open-time ex- 
tragalactic key project on the observatory, and will eventu- 
ally observe ~550 sq deg in five p assbands at 100, 160 , 250, 
350 and 500 p with the PACS (|Poglitsch et alJl201Ch and 
the SPIRE (|Griffin et al.1 12010| ) instruments, allowing the 
study of interstellar dust at angular scales ~ 10 times smaller 
than with IRAS. With these data, approximately 250,000 
galaxies out to redshifts of 3 to 4 are expect ed to be detected. 
We omit the PACS data jlbar et alj|2010h from our analy- 
sis as the selected sources may possibly be contaminated 
by warmer dust heated by the AGN and the data do not 
reach the required depth to be useful. Instead, we h ave con- 
centrated on the SPIRE data (|Pascale et alj|201ll ). SPIRE 
wavelengths have been chosen because they are less likely to 
be contaminated by direct emission from the warm dust of 
the AGN torus, with a starburst component at these wave- 
lengths being required to acco unt for the total FIR emission , 
as sugges ted by the SEP fits oflHatziminaoglou et al.l (|2010T l . 

The iGibson et all (|2009h catalogue has been cross- 
matched with the H-ATLAS 9, 12 and 15 h fields. Each field 
is approximately 12° in RA by 3° in Dec (6° by 3° for the 
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Figure 1. Distributions of BAL and non-BAL QSOs as a func- 
tion of redshift and absolute magnitude in SDSS i band. The top 
panel shows the full sample before our selection methods have 
been used. The dashed box shows the redshift and absolute mag- 
nitude cut-offs used in this study. The lower panel shows the final 
sample after removal of those QSOs falling outside our absolute 
magnitude and redshift cutoffs, and after statistical matching (see 
text for details). Red asterisks are non-BAL QSOs, blue diamonds 
are HiBAL QSOs, the black filled star is our solitary LoBAL QSO, 
green circles are our detected HiBAL QSOs and orange squares 
are our detected non-BAL QSOs. 



other properties of our sample. As a result of this distinction, 
the total number of BAL QSOs following the 'classic' BI def- 
inition drops to 36. One of these, SDSS115404.13+001419.6, 
is a LoBAL QSO. This source is not considered in our stack- 
ing analyses in Section 4. The distribution of BAL QSOs 
over the 3 fields is as follows: 9 h field - 19, 12 h field - 12, 15 h 
field - 19. Our comparison sample consists of 329 non-BAL 
QSOs drawn from the SDSS Data Release 5 lSchneider et all 
(|2007f ). We did not use the SDSS DR8 catalogue since there 
has been no corresponding BAL QSO catalogue produced. 
We trimmed objects at random from overpopulated parts of 
the comparison sample (determined via visual inspection of 
Mi and z) until a two-sample KS test comparison in Mi and 
z between the non-BAL and BAL QSO sample returned a 
null-hypothesis probablility of p > 0.05; for the final non- 
BAL QSO comparison sample separate KS tests on redshift 
and Mi give p- values of 0.67 and 0.20 respectively ('extended 
sample') and 0.89 and 0.17 respective ly ('classic' sample). 
Running a 2-d KS test (|Peacock|[l983l ) on redshift and Mi 
for both samples returns a p-value of 0.27. We can therefore 
assume the populations are adequately matched in Mi and 
z. The lower panel of Fig. [1] shows the final BAL and non- 
BAL QSO samples on the (z, Mi) plane while Table [T]gives 
details of the final BAL QSO sample. 

In the next section, we compare the FIR properties 
of matched BAL and non-BAL QSO samples, and deter- 
mine SFRs for the samples based on separate classification 
schemes (as explained in Section 2). Our study improves on 
previous work since whilst the H- ATLAS survey is still not 
complete, the three fields in the H-ATLAS phase 1 dataset 
centred at 9, 12, and 15h give us a large area providing a 
uniformly selected sample of objects cross-matched with the 
Sloan Digital Sky Survey (SDSS). 



4 METHOD AND RESULTS 

Here we determine the FIR properties of the selected QSOs 
and investigate how the BAL QSOs relate to the non-BAL 
QSOs. We further discuss the effects of different selection 
criteria i.e., BI versus BIq. 



12 h field), which places a limit on the size of our final sam- 
ple. Only those QSOs falling within the field boundaries are 
included resulting in an initial sample of 88. These further 
reduce to 83 due to 5 that have a C IV Balnicity Index value 
of 0. These BAL QSOs have been classified using alternative 
lines such as Siiv or Mgn, and are considered beyond the 
scope of this paper. 

We next apply constraints within the parameter space 
1.5 < z < 2.3 and -28 M sC -24, where most of the 
BAL QSOs lie, in order to reduce the susceptibility of our 
study to the effects of luminosity or redshift evolution (e.g. 
iBonfield et all 1 20 111 ). This restriction gives two advantages: 
it allows easy construction of a non-BAL QSO comparison 
sample and limits the effects of cosmological evolution whilst 
comparing similar luminosities. In this way the sample is 
further reduced to 50 BAL QSOs. We have also addressed 
the important distinction between the 'classic ' BI and the 
'extended' version given in iGibson et all |2009l ). Bio, to see 
if this classification has any effect on the flux densities and 



4.1 Stacking analysis 

The FIR flux densities used in this section are taken directly 
from the final PSF-convolved images for all three H-ATLAS 
fields. Cutouts of set size (210x210 arcsec 2 ) of the region 
surrounding each BAL/non-BAL QSO are generated, and 
these are co-added, each pixel being the weighted mean flux 
density in that pixel. The weighted mean stacked flux den- 
sity value of each source is taken to be that at the centre 
of the stacked image which i s the closest p i xel to the cat- 
alogue position in eit her the IGibson et al.l l|2009f l BAL or 
ISchneider et~ai] (|2007l ) SDSS QSO samples. We background 
subtract by randomly selecting a sample of 40,000 pixels 
from each field to create a mean background for that field 
at that wavelength. We then subtract this mean background 
from every pixel in our cutout images for each QSO. These 
mean values for each wavelength are shown in Table [2] while 
Fig.[2]shows the BAL and non-BAL QSO stacked images at 
each wavelength. 

We find that at longer wavelengths the stacked BAL 
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Table 1. The BAL QSO Sample. The classification method used to determine the Balnicity of a source is given in the Type/Detected 
column, along with whether it is detected. If a source is classified using BI, then it will also be classified as a BAL QSO by the BIo 
definition. The 250 fim flux density and signal-to-noise ratio are given derived from the 250 Ami fl ux de nsity divided by the sum in 
quadrature of the instrumental noise in the corresponding pixel and confusion noise in iRigbv et al, ] J201l|y 5<r detections found in the 
Phase 1 catalogue with matched optical/NIR counterparts (see Section 4.2) are shown with their flux density values taken from the 
Phase f catalogue. The average 5<r limits at 250, 350 and 500 (im are 33.5, 37.7 and 44.0 mjy beam -1 . For undetected sources, their 
flux density values are the best estimates from measuring the flux density in the closest pixel of the PSF-smoothed map and noise map 
and are not reliably associated with R > 0.80 to a 250 fim source as for our detections. The FIR data for the non-BAL Quasars within 
H-ATLAS will be dealt with in a future paper (Bonfield et al. in prep), however all maps and catalogues used in this paper will be 
released as part of the Phase 1 data release from the H-ATLAS website. 



Source z Type/ Civ EW M; 250 fim flux SNR 350 fim flux SNR 500 fim flux SNR 

Detected? (A) density (mjy) 250 density (mjy) 350 density (mjy) 500 
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QSO image is less distinguishable from the background, 
which is expected due to increased confusion noise and a 
higher instrumental noise; galaxies at 500 fim are also ex- 
pected to be fainter, since the Rayleigh- Jeans tail is being 
sampled. Table[2]shows the central pixel values of the stacks. 
We find that there is no statistical difference between BAL 
QSOs and non-BAL QSOs; their flux densities in all band- 
passes are the same within the errors. This result holds for 
both the 'classic' and 'extended' BAL QSO samples. 

We attempt to reproduce our results using the Institut 
d'Astrophysique Spatiale (IAS) stacking library. It stacks 



data to allow a statistical detection of a faint signal using 
positions of galaxies detected at shorter wavelengths as with 
our method but has also been tested and validated with 
galaxy clus tering in mind, something which may affect our 
results (see lBavouze t 2008; B cthermin et al.ll2010l for further 
details). However, we still find no difference between the 
BAL QSO and non-BAL QSO samples in each bandpass 
and the IAS values are in agreement with our values within 
the errors (see Table |2J. 

It should be noted that submillimetre fluxes may still 
be overestimated due to clustering of sources if they emit in 
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250 |im 350 fim 500 fan 

Non-BAL Non-BAL Non-BAL 

329 sources 329 sources 329 sources 



Figure 2. The 250, 350 and 500 /an weighted mean flux density stacks of the 'extended' BAL and non-BAL QSOs (top and bottom row 
respectively) samples. Each postage stamp image is 210x210 arcsec 2 . In each corner in red is given the wavelength, classification, and 
the number of sources used to create that stack. 



the SPIRE bands, but since we cannot identify the sources 
around the QSO at high redshifts there is little that can be 
done to solve this problem. We also emphasize that the ef- 
fects of confusion should be the same for both BAL QSOs 
and non-BAL QSOs since they are treated in an equal man- 
ner within both methods. 

The Gaussian errors quoted in Table [2] must be treated 
with caution, since background noise will be non-Gaussian 
owing to confusion noise. To determine quantitatively 
whether the stacked BAL QSOs are detected significantly, 
flux densities from 120,000 random ly chosen positions in th e 
field have been measured following lHardcastle et all (|201Cf ) . 
The random background compared with the BAL and non- 
BAL QSO fluxes are shown in Figure [3] Using a KS test, 
we can examine whether the flux densities of the stacks are 
statistically distinguishable from those taken from randomly 
chosen positions, as a KS test is not influenced by the noise 
properties. 

In comparison to the background we detect the BAL 
QSOs at all wavelengths using both classification schemes. 
The significance of the 250 )im detection for 'extended' BIo 
BAL QSOs is well over 5a (p < 10~ 10 ), and as expected, 
the significance decreases towards longer wavelengths, being 
lowest at 500 /im (p~5xl0~ 4 ). The larger sample of non- 
BAL QSOs is also detected with high significance above the 
background at all wavelengths (p < 10~ 10 ). Returned KS 
probabilities are shown in Table[3] We also used a KS test to 
compare the distribution of flux densities between the BAL 
QSO samples and the non-BAL QSO samples (see Fig. [3j. 
At 250, 350 and 500 their p- values are p = 0.39, 0.80, and 
0.37 respectively for the 'extended' sample when comparing 
the BAL QSO sample with the full non-BAL QSO sample. 
We therefore cannot reject the null hypothesis that the BAL 
and non-BAL QSOs are drawn from the same underlying 
flux density distibution. 

We performed a second test where the 120,000 stacks 
were created from flux densities extracted from random po- 
sitions, with the number of elements in the stack equal to 
the elements in the comparison sample stack, i.e. 49, 35 and 



329, and we then compared these random stacks with the 
stacked flux densities of our BAL and non-BAL QSO sam- 
ples. The fraction of random stacks where the weighted mean 
flux density exceeds the weighted mean flux density in our 
sample stacks provides an estimate of the probability of the 
detection of our stacks being just a background fluctuation. 
We find that 250 fim is generally of greater significance than 
that at 500 ^m for the BAL QSO samples (e.g. p < 10" 5 , 
p = 4.8 x 10~ 4 respectively for the 'extended' BI BAL QSOs 
in each bandpass), and that none of the 120,000 flux density 
stacks extracted at random positions have larger weighted 
mean flux densities than the BAL and non-BAL stacked 
flux density as shown in Table [3] This indicates that our 
weighted mean flux density values in Table [2] are not due to 
background fluctuations. 



4.2 Detection rates 

We now investigate how many individual BAL QSOs are 
significantly detected in the H- ATLAS data at the >5<r 
level. The Phase 1 H-ATLAS catalogue used a likelihood 
ratio technique to identify optical and near-infrared coun- 
terparts wi th objects detecte d by Herschel as with the SDP 
catalogue (ISmith et al.l l201lh in the phase 1 data. We find 
that three BAL QSOs from our sample are associated with 
(5 a) 250 fim sources in the H- ATLAS Phase 1 catalogue (i.e. 
they have reliability, R > 0.80); two of these are "classic" 
BJ-defined BAL QSOs. For the non-BAL QSOs, 27 of the 
329 have reliable R > 0.80 250 /im counterparts. This raises 
the question of whether there is an actual difference in the 
detection rate of BAL QSOs and non-BAL QSOs. We de- 
tect 6.0 per cent of the 'extended' sample (8.3 per cent of the 
'classic' BI defined sample) compared to 8.2 per cent of the 
non-BAL sample. Comparing these detection rates using a 
binomial probability distribution gives a null result (p — 0.19 
and p = 0.23 for the 'extended' and 'classic' sample respec- 
tively). We cannot therefore reject the null hypothesis that 
the detected HiBALs are drawn from the same distribution 
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Table 2. The BAL ('extended' and 'classic' samples) and non-BAL QSO FIR weighted mean flux densities in the 250, 350 and 500 /im 
bandpasses. The number in brackets gives the number of objects within each stack. We also include here our determined background 
noise for each field at each wavelength. 
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Table 3. The KS statistics and probabilities of each sample being indistinguishable from a randomly selected sample of flux densities 
taken from the H-ATLAS fields. The fraction shows how many random stacks had flux densities greater than our sample stacks. 
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as the detected non-BAL QSOs, consistent with our conclu- 
sions above on their flux density distribution. 

4.3 FIR luminosities and star-formation rates 

We have limited information on the shape of the SED of 
our QSOs and therefore must choose a suitable template 
in order to compute FIR luminosities (Lfir). The primary 
galaxy that we have chosen for this purpose is Mrk 231. Of 
the local ULIRGs, it is the most luminous, wit h a rest-frame 
lumin osity at 8 - 1000 fj,m of 3.2 x 10 12 L (|Sanders et al.l 
I2003T I, and it shows similar broad absorption line features to 
our sample, havi ng been classified as a BAL QSO, strictly 
a LoBAL QSO (|Smith et al.lll995l ). We have also investi- 
gated the effect of using a template based on IZwl, a local 
Seyfert type 1 ga laxy in the nearby universe at z — 0.0589 
|Ho fc Kirnll2009l ). We assume T = 45 K and emissivity in- 
dex P = 1.6 for Mrk 231, and T = 36 K and p = 1.7 for 
IZwl from fitting a greybody to data points retrieved from 
the NASA Extragalactic Database (32 and 9 respectively). 

We determine the FIR luminosity via a ratio method 
between the 250 /im flux density of our greybody template 
if placed at the redshift of the QSO in our sample and the 
250 /im flux density of the QSO. This implicitly assumes 
that the QSO has the same FIR SED as the greybody tem- 
plate. The calculation of the SFR was then performed in the 
st andard manner us ing the relationship with Lfir published 
by iKennicutt] (|l99&f ) in equation [2] and these are presented 
in Table H i.e. 

SFRiMeyr' 1 ) = 4.5 x 10~ 44 L F /fl (ergs- 1 ). (2) 

ThelKennicuttl Jl998l) relatio n is dependent on the models of 
iLeitherer fc Heckmanl ll 19951 ) . These models were specifically 
created to correspond to the conditions that are prevalent 
within giant Hn regions, Hn galaxies, nuclear starbursts and 
infrared luminous starbu rst galaxie s and traces diffuse gas 
within these regions. The lKennicuttl (| 19981 ) relation assumes 
continuous starbursts of age 10 — 100 Myr, and requires the 
integrated IR luminosity over the range 8 — 1000 ^m. 

The average Kennicutt SFR is 240 ± 21 M yr" 1 for 
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Flux in Pixel (mjy) 

Figure 3. The distribution of flux density at each wavelength for 
the full BAL (blue outline) and non-BAL QSO (green outline) 
samples in bins of width 1 mjy. The random background flux 
densities are shown as a red outline. 



the 'extended' and 223 ± 25 Mq yr~ for the 'classic' sam- 
ple (Table [4); these values are consistent. Note that these 
estimates are determined from weighted stacking of the in- 
dividual luminosities and SFRs derived from the 250 /im 
flux density of each galaxy rather than calculated from 
the weighted stacked flux density values in Table [2] The 
BAL and non-BAL average SFRs are also consistent within 
the errors. We find that the 3 detected BAL QSOs have 
SFR > 1000 Mqvt - . Since we cannot exclude the possi- 
bility that some fraction of the FIR emission comes from 
dust heated by the AGN, these SFRs should be regarded as 
upper limits. We also find that the average FIR luminosities 
presented in Table 0] for both the BAL and non-BAL QSO 
samples are, on their own, sufficiently large to classifiy them 
as ULIRGs. We note that the derived SFRs for each SPIRE 
band for IZwl are more consistent than those determined for 
Mrk 231. This suggests that dust in BAL QSOs at higher z 
is likely cooler than that found in the low-z analogue Mrk 
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Table 4. The star-formation rates and FIR luminosities of the detected (marked with an asterisk) BAL QSOs within the sample, along 
with the average star-formation rate in the 250, 350 and 500 /an bandpasses and weighted mean for the 'extended' , 'classic' and non- 
BAL samples as a whole using Mrk 231 and IZwl as templates. We also include here our rare LoBAL QSO 115404.13+001419.6. We 
emphasize however that it is not reliably associated with R > 0.80 to a 250 (im source as with our detections and instead is the best 
estimate from measuring the flux density in the closest pixel of the PSF-smoothed map. At 500 /im it has a negative flux so SFRs are 
not included. The quoted Lfir (using our ratio method for 250 fim) and SFR estimates are determined from weighted stacking of the 
individual luminosities and SFRs of each galaxy rather than calculated from the weighted stack values in Table [2] (shown in the flux 
density column), though the error will be underestimated as a result of confusion noise. 
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Figure 4. The di stribution of C IV a bsorption-line equivalent 
width for the full I Gibson et al.l | |2009|) sample (grey) and our 
'classic' BAL QSO sample (blac k outli ne) in bins of width 5 A. 
One source in the IGibson et ail ||2009|) sample is omitted here, 
SDSS090901. 80+535935.7, having a Civ absorption-line equiva- 
lent width of 193.9 A. All others have Civ absorption-line equiv- 
alent widths below 90 A. 



231. Finally, we have also investigated the effect of varying 
the dust temperature and j3 values and find that for plau- 
sible values (25 K < T < 60 K and 1.1 < /3 < 1.8) there is 
about a factor of 2 variation in the final derived SFR values 
which does not change our conclusions. 



4.4 FIR luminosity and C IV absorption— line 
equivalent width 

From an analysis of submillimetre wavelength data, 
iPriddev et ail (|2007l ) found a weak correlation between FIR 
flux density and C IV absorption-line equivalent width and a 
rather stronger statistical link between submillimetre detec- 
tion rate and Civ absorption-line equivalent width; i.e., of 
the 15 sources they observed, all 6 (>2<r) detections in the 
850 fim band had a C iv absorption-line equivalent width of 
> 25 A. A KS test returned a probability of 0.01 that the de- 



tections and non-detections had the same C iv absorption- 
line equivalent width. 

The small number of detections in our sample means we 
cannot perform the same tests that Priddey et al. used. In 
our sample there are only 5 sources with C iv absorption- 
line equivalent widths of > 25 A in the 250 /an bandpass, 
of which 3 are detections. Assuming a binomial probability 
distribution we estimate the probability of detections hav- 
ing C IV absorption-line equivalent widths > 25 A by chance 
is only 0.005. However, these are detections in one band- 
pass alone and the number of objects with C iv absorption- 
line equivalent width > 25 A is very small. We must there- 
fore choose a different C iv absorption-line equivalent width 
cut-off for our HiBAL QSOs, such that we have similar 
sized samples to compare. We set this cut-off at ^ 20 A, 
which gives us 11 HiBAL QSOs to compare with the re- 
maining 24 in the classic sample. For completeness we tested 
whether the distribution of C iv absorption-line equivalent 
widths within our restricted BAL QSO sample was represen- 
tative of thos e with in the parent population as described by 
IGibson et all (|2009h . A KS test on the C iv absorption-line 
equivalent widths of each population returns p = 0.12, so 
we are able to assume our sample is representative. Figure^] 
shows histograms of the two C iv absorption-line equivalent 
width distributions. For this analysis we studied the 'clas- 
sic' BAL QSO sample only since it suffers from less con- 
tamination from mis-classified sources than the 'extended' 
sample due to the strictness of the BI mechanism. Further- 
more, any true underlying physical relationship involving 
C iv absorption-line equivalent width will be expected to 
correlate with luminosity rather than flux density. We there- 
fore determined the Lfir luminosities of each BAL QSO 
using our assumed T and j3 for Mrk 231 and IZwl, sep- 
arated them using o ur 20 A cutoff and then performed a 
Mann- Whitney test (|Mann fc Whitnevl[l947l ) to determine 
whether a greater luminosity correlates with a larger C iv 
absorption-line equivalent width. For the Mrk 231 and IZwl 
templates respectively, the test returned probabilities of 0.23 
and 0.22, which indicates that there is no significant evidence 
for a difference in the two subsamples. Obviously with this 
method there is the caveat that the BAL QSOs must each 
have a similar SED shape to our greybody templates but 



H- ATLAS: The FIR Properties of BAL Quasars 9 



these results indicate there is no evidence for a dependence 
of FIR luminosity on Civ EW. 



5 DISCUSSION 

Our study shows that the FIR properties of HiBAL QSOs 
are statistically indistinguishable from those of non-BAL 
QSOs. We also find that FIR emission has no dependence 
on C IV absorption-line equivalent width. Therefore the FIR 
luminosity has no dependence on the absorption strength of 
the outflow, and our results do not require an evolutionary 
link between HiBAL and non-BAL QSOs. Rather a simple 
orie ntation effect argume nt is sufficient. 

iPriddev et aL I (|2007f ) found tentative evidence for a link 
between C IV absorption-line equivalent width and submil- 
limetre wavelength detection (albeit with a 2a detection 
threshold) . As discussed by Priddey et al. such a link could 
be used to argue that BAL QSOs are i n a distinct evolu- 
tiona ry phase along the lines proposed bv lPage et al.l (|2004l . 
l201lf ) where the QSO wind is in the process of terminat- 
ing an epoch of enhanced star formation. Our results do not 
support this weak evidence for such a link between FIR out- 
put and C IV absorption-line equivalent width. The discep- 
ancy may be linked to a possible selection effect discussed 
by Priddey et al.; that the most reddened QSOs must be 
intrinsically more luminous to meet their blue selection cri- 
terion, and therefore the most extreme objects, which might 
host the most powerful outflows would have higher dust red- 
dening and higher FIR/submm output. This could explain 
why we see no such correlation, since the SDSS i band does 
not suffer reddening to as large a degree as the bluer bands 
in which past samples of QSOs were selected. 

However, a caveat to our results is that our sample is 
composed almost entirely of HiBAL QSOs. We have not con- 
sidered either LoBAL QSOs or their rarer FeLoBAL cousins, 
which may well show FIR properties inconsistent with a 
simple orientation s cheme hypothesis (s ee be low). We note 
that t he samples of IWillott et al.l (|2003h and IPriddev et all 
were a l so m ostly composed of HiBAL QSOs. The 
IPriddev et al.l (|2007t ) s ample contain s only 1 known LoBAL 
QSO, and many of the lWillott et all (|2003l ) QSOs classified 
as LoBAL QSOs are uncertain due to Mg II moving out of 
the spectral range at z > 2.26. 

Regarding studie s at other wavelengths, the findings 
of [G rccn et all (|200ll ) carried out using the Chandra X- 
ray observatory provide weak evidence that HiBAL QSOs 
appear X-ray weak due to intrinsic absorption, and that 
their underlying emission is actually the same as that of 
non-B AL QSOs (|Gallagher et alj|2002t l200d : iGibson et ail 
120091 ), which supports an orientation scheme. In contrast, 
even taking into account intrinsic absorption, LoBAL QSOs 
remain X- ray weak. An anal y sis of X-ray absorption per- 
formed bv lStreblvanska et all (|201Ch suggests that LoBALs 
and HiBALs may be physically different objects. This is re- 
inforced by c omposite SEDs r angin g from X-ray to radio 
generated by iGallagher et all (|2007l ) who find similarities 
between the SEDs of HiBALs and non-BAL QSO but dif- 
ferences between these populations and LoBALs. However, 
it should be noted that the absorbed nature of BAL QSOs 
means that their X-ray data are limited in quality, so these 
findings cannot be taken as conclusive. 



Certainly LoBAL QSOs are found to have redder spec- 
tra than other types |Reichard et all l2003bl ) , implying a 
larger dust mass in the host galaxy and/or vicinity of the 
SMBH. The idea that LoBALs and FeLoBALs form a sep- 
arate evolutionary class finds suppor t from the work of 
IFarrah et all (|2007T ) and lllrrutia et all (|2009h . IFarrah et all 
|2007l ) argue that FeLoBAL QSOs are galaxies where a mas- 
sive merger driven starburst is ending, and the last remnants 
of a dust cocoon are being removed by a rapidly accreting 
supermassive black hole at the galaxy's centre. The SED fits 
of those detected at longer wavelengths (4/9) require a star- 
burst component of the order of several hundred solar masses 
per year, which is taken to imply that FeLoBAL QSOs are 
associated w i th UL IRGs. Similar conclusions are drawn in 
IFarrah et all (|2010t ) who find that FeLoBAL QSOs span a 
range of spectral shapes, consistent with the idea that they 
may be a transition phase in the lifetime of an AGN. How- 
ever, as noted in their papers, the small and inhomogeneous 
sample sizes render s uch conclusions tent ative. A larger sam- 
ple was presented by IFarrah et all | |2012] ) of 31 objects with 
optical to far-infrared photometry. These were found to all 
be highly luminous (> 10 12 L©), yet the bulk of IR emission 
in the majority of sources came from the AGN rather than 
a starburst component. The mid-IR and far-infrared prop- 
erties of the FeLoBAL QSOs spanning the redshift range 
0.8 < z < 1.8 presented by Farrah et al. further reinforce 
the idea that they are a class apart from the general QSO 
pop ulation. 

lUrrutia et all i|2009l ) also found an unusually high per- 
centage (32 per cent using the conservative BI selection) of 
LoBAL QSOs in a sample of dust reddened type-1 QSOs; 
in fact all of these objects can also be classified as FeLoB- 
ALs. The orientation hypothesis is irreconcilable with such 
findings allowing an evolutionar y hypothesis to be in voked, 
even with the caveats offe red by Urrutia et all (|2009r ) . 

The subset of the full lGibson et alH 20091 ) sample stud- 
ied here contains only one QSO, SDSS115404.13 + 001419.6, 
classified as a LoBAL. These objects are rare, making 
up only 15 per cent of all observed BAL QSOs. A very 
deep C IV absorption trough (greater than anything seen 
in the HiBAL sample), combined with a high flux density 
(29.88 ± 6.65 mjy at 250 (im) makes it a conspicous member 
of our sample. However, little can be said of an entire popula- 
tion from one object. Future work in the FIR/submillimetre 
wavebands targeting LoBALs and/or FeLoBALs would be 
invaluable. 



6 CONCLUSIONS 

(i) Using a stacking analysis we have determined that 
HiBAL QSOs at 1.5 ^ z < 2.3 are statistically indistin- 
guishable in terms of FIR luminosity (and therefore SFR) 
from a matched sample of non-BAL QSOs. This result is 
broadly consistent with previous work conducted at sub- 
millim etre wavelengths (|Willott et all 120031 ; IPriddev et all 
The average FIR luminosities of both our HiBAL and 
non-BAL samples are > 10 12 Lq, sufficient to classify them 
as ULIRGs. We calculate weighted mean SFRs (strictly up- 
per limits using this relation due to the possibility of dust 
heating by the AGN and our BAL QSOs having overdensi- 
ties of neighbours) for the HiBALs of 240 ± 21 M Q yr" 1 for 
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the 'extended' and 223 ±25 Mq yr 1 for the 'classic' sample 
from the flu x density at 250 fim. 

(ii) While IPriddev et all <|2007h found tentative evidence 
for a dependence of submillimetre flux density on C IV 
absorption-line equivalent width, we find no such depen- 
dence at FIR wavelengths. We suggest that the Priddey et 
al. result may have been due to small number statistics and 
selection effects as noted in the discussion. 

(iii) Within our samples, 3/49 HiBALs and 27/329 non- 
BALs are detected at > 5 a significance. The detection rates 
for the two species are statistically indistinguishable. 

(iv) Taken together these results suggest that BALs 
(strictly HiBALs) can be unified with non-BAL QSOs within 
a simple orientation scheme where a BAL QSO is observed 
only if a nuclear outflow intercepts our line-of-sight. How- 
ever, with the current data we are unable to say whether 
(Fe)LoBALs can similarly be accommodated within this 
scheme or whether they form a distinct population perhaps 
caught at a key phase in their evolution. Future observations 
at FIR/submillimetre wavelengths would be valuable in this 
respect. 
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